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1 Introduction 

"It seems there is no problem in modern physics for which there are on record 
as many false starts, and as many theories which overlook some essential fea- 
ture, as in the problem of the thermal conductivity of nonconducting crystals^^ 
[1]. This statement by R. Peierls goes back to almost 50 years ago, yet it ap- 
pears to be still valid. Compared with charge flow (electric current), much 
less is known about the heat flow. 

Nevertheless, encouraging results have been obtained in recent years (see 
reviews [2] and the references therein). For example it is now known that in 
one dimensional systems of the Fermi-Pasta- Ulam (FPU) type [3], heat con- 
duction is anomalous and the coefficient of thermal conductivity k diverges 
with the system size L as k ~ L^^^ (when the transverse motion is con- 
sidered K ~ L^^^ [4]). The connection between anomalous conductivity and 
anomalous diffusion has been also established[5, 6], which implies in partic- 
ular that a subdiffusive system is an insulator in the thermodynamic limit 
and a ballistic system is a perfect thermal conductor, the Fourier law be- 
ing therefore valid only when phonons undergo a normal diffusive motion. 
More profoundly, it has been clarified that exponential dynamical instability 
is a sufficient [7, 8] but not a necessary condition for the validity of Fourier 
law [9, 10, 11, 12]. These basic studies not only enrich our knowledge of the 
fundamental transport laws in statistical mechanics, but also open the way 
for applications such as designing novel thermal materials and/or devices 
such as the thermal rectifier[13, 14] and the thermal transistor [15]. 

In this paper we give a brief review of the relation between microscopic 
dynamical properties and the Fourier law of heat conduction as well as the 
connection between anomalous conduction and anomalous diffusion. We then 
discuss the possibility to control the heat flow. 
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2 Microscopic dynamics and the Fourier heat law 

In spite of several years of investigation, so far there is neither phenomeno- 
logical nor fundamental transport theory which can predict whether or not a 
given classical, many-body Hamiltonian system, yields an energy transport 
governed by the Fourier heat law. It is known that heat flow obeys a simple 
diffusion equation which can be regarded as the continuum limit of a discrete 
random walk. However, modern ergodic theory tells us that for K-systems, 
a sequence of measurements with finite precision mimics a truly random se- 
quence and therefore these systems appear precisely those deterministically 
random systems tacitly required by transport theory. It is therefore interest- 
ing to establish which class, if any, of many-body systems satisfy the necessary 
stringent requirements for the validity of Fourier law of heat conduction. 

2.1 Dynamical chaos is a sufHcient condition for the Fourier law 
Ding-a-ling model 

The first example for which convincing evidence has been provided that 
Fourier law can be derived on purely dynamical grounds, without any ad- 
ditional statistical assumptions, is the so-called ding-a-ling model proposed 
in [7]. 

The Hamiltonian of the ding-a-ling model is: 

^^\Y.(Pk+ '^I'il) + ^ard point core (1) 

where Uk equals to for even k and zero for odd k and where all particles 
have unit mass. It can be shown that it is possible to fix the energy per 
particle {E — 1) and the half of the lattice distance between two bound 
particles (/q = 1). After this the dynamics is uniquely determined by the 
frequency u> which becomes therefore the only adjustable parameter. 



1 2 3 4 N-1 N 

Fig. 1. The N-particle ding-a-ling model. Odd particles can move freely in between 
two collisions, while even particles are bounded by a harmonic spring. 



As it seen in Fig 1. the model is a one-dimensional array of equal mass, 
hard-point particles, the even-numbered particles form a set of equally spaced 
harmonic oscillators with each oscillator bound to its individual lattice site 
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and with all oscillators vibrating with the same frequency oj. The odd- 
numbcrcd particles are free particles constrained only by the two adjacent 
even- numbered oscillators. 

This model has been chosen in order to meet two requirements. First, 
it is necessary to select a dctcrministically random system, and second, the 
system must be sufficiently simple to allow efficient numerical analysis. Let 
us also recall that even systems obeying the Fourier heat law can transport 
energy in the form of slowly decaying coherent excitations such as soundlikc 
pulses. In numerical experiments, which unavoidably consider only a small 
number of particles, this phenomenon is quite troublesome. Though K or 
almost K-systcm guarantees that these soundlikc solutions will eventually 
decay, one needs to find a small chaotic system in which this decay rate is 
sufficiently rapid[16]. The selected model (1) is a many-body system which 
exhibits, as the frequency is varied, the full range of behaviour from intcgrable 
to almost K and which at the same time has no problem with energy-bearing, 
long-lived, solitonlike pulses. Indeed, it can be shown that as u) is increased 
from zero (which is the well known intcgrable 1-d hard point gas), the system 
undergoes a transition from intcgrable to almost fully chaotic. This fact makes 
this model ideal to study the connection between Fourier law and microscopic 
dynamical chaos. 

Heat conductivity has been studied by placing the end particles in contact 
with two thermal reservoirs at different temperatures (see [7] for details)and 
then integrating the equations of motion. Numerical results [7] demonstrated 
that, in the small u regime, the heat conductivity is system size dependent, 
while at large w, when the system becomes almost fully chaotic, the heat 
conductivity becomes independent of the system size (if the size is large 
enough). This means that Fourier law is obeyed in the chaotic regime. 

Lorentz gas channel 

The above conclusions have been nicely confirmed and clarified in [8] where 
the heat conduction has been studied in a Lorentz gas channel - a quasi 
one dimensional billiard with circular scatterers. The model (see Fig. 2) con- 
sists of two parallel lines and a series of semicircles of radius R placed in a 
triangular lattice along the channel. By construction no particle can move 
along the horizontal direction without colliding with the disks. The dynam- 
ics in the Lorentz gas is rigorously known to be mixing and all trajectories 
with nonzero projection on the horizontal direction are of hyperbolic type; 
further it has positive Kolmogorov-Sinai entropy and a well defined diffu- 
sive behavior. Very accurate numerical evidence has been provided [8] which 
shows that heat conduction in this model obeys the Fourier law. This result 
therefore clearly indicates that mixing with positive Lyapounov exponents is 
a sufficient condition to ensure Fourier law of heat conduction. 
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Fig. 2. Tlie geometry of tiie Lorentz gas ciiannel model. Tlie two lieat reservoirs 
at temperatures Tl and Tr are indicated. 



2.2 Chaos is not a necessary condition 
Triangle billiard gas channel 

Quite naturally, the next question which arises is whether strong, exponential 
unstable chaos, being sufficient, is also necessary. 

In this connection let us remark that in spite of several efforts, the rela- 
tion between Lyapounov exponents, correlations decay, diffusive and trans- 
port properties is still not completely clear. For example a model has been 
presented [17] which has zero Lyapounov exponent and yet it exhibits un- 
bounded Gaussian diffusive behavior. Since diffusive behavior is at the root 
of normal heat transport then the above result [17] constitutes a strong sug- 
gestion that normal heat conduction can take place even without the strong 
requirement of exponential instability. 




Fig. 3. The geometry of the triangle billiard gas channel. Particles move in the 
region outside the triangular scatterers. The two heat reservoirs at temperatures 
Tl and Tr are indicated. 



The above problem has been addressed in [11], where we have considered 
a quasi-one dimensional billiard model which consists of two parallel lines 
and a series of triangular scatterers (see Fig. 3). In this geometry, no parti- 
cle can move between the two reservoirs without suffering elastic collisions 
with the triangles. Therefore this model is analogous to the previous Lorentz 
gas channel with triangles instead of discs: the essential difference is that 
in the triangular model the dynamical instability is linear and therefore the 
Lyapounov exponent is zero. 

Strong numerical evidence has been given [18] that the motion inside 
a triangular billiard, with all angles irrational with tt is mixing, without 
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any time scale. It is therefore reasonable to expect that the motion inside 
the irrational polygonal area of Fig 3 is diffusive thus leading to normal 
conductivity. 

Indeed, numerical results in [11] show that in the irrational case (when the 
ratio O/n and ^/tt are irrational numbers) the system in Fig 3 exhibits normal 
diffusion and the heat conduction obeys the Fourier law. In the rational case 
instead, the system shows a superdiffusive behavior, (cr^) = 2Dt^'^'^^ [lljand 
the heat conductivity diverges with the system size as k ~ ^o.25±o.oi^ 

Alternate mass-core hard potential channel 

In the two billiard gas models just discussed there is no local thermal equi- 
librium. Even though the internal temperature can be clearly defined at any 
position[8], the above property may be considered unsatisfactory [19]. In or- 
der to overcome this problem, we have recently introduced a similar model 
which however exhibits local thermal equilibrium, normal diffusion, and zero 
Lyapunov exponent [12]. 
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Fig. 4. The geometry of the alternate mass-core hard potential channel. The ele- 
mentary cell is indicated by the two dotted lines. The bars have mass M = 1, and 
the particles have mass m = {\/5 — l)/2. The two heat baths at temperatures Tl 
and Tr are indicated. 



This model consists of a one-dimensional chain of elastically colliding par- 
ticles with alternate masses m and M . In order to prevent total momentum 
conservation we confine the motion of particles of mass M (bars) inside sep- 
arate cells. Schematically the model is shown in Fig. 4: particles with mass m 
move horizontally and collide with bars of mass M which, besides suffering 
collisions with the particles, are elastically reflected back at the edges of their 
cells. In between collisions, particles and bars move freely. 

Numerical results [12] clearly indicate that this model also obeys the 
Fourier law. 

3 Anomalous Heat Conduction 

Numerical experiments have shown that in many one dimensional systems 
with total momentum conservation, the heat conduction does not obey the 
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Fourier law and the heat conductivity depends on the system size. For exam- 
ple, in the so-called FPU model, k ^ L^, with /3 = 2/5, and if the transverse 
motion is introduced, /3 = 1/3. Moreover, in the billiard gas channels (with 
conserved total momentum), the value of P differs from model to model [5]. 
The question is whether one can relate (3 to the dynamical and statistical 
properties of the system. 



Fig. 5. Comparison of prediction (4) with numerical data. Normal diffusion (★). The 
ballistic motion (*). Superdiffusion: ID Ehrenfest gas channel[9](v); the rational 
triangle channel[ll] (empty box); the polygonal billiard channel with 4>i = (v^ — 
l)7r/4), and 02 ~ 7r/3 [10] (A); the triangle-square channel gas[6] (o). /3 values 
are obtained from system size L £ [192, 384] for all channels except Ehrenfest 
channel[9]. The FPU lattice model at high temperature regime[6](0), and the single 
walled nanotubes at room temperature (®). Subdiffusion: model from Ref.[10] (solid 
left triangle). The solid curve is l3 = 2 — 2/ a. 



Recently, a simple formula has been found [5] which connects anomalous 
heat conductivity with anomalous diffusion. More precisely, it has been shown 
that for a one dimensional system, if the energy diffusion can be described 




l.Or 
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a 



by 



(0-2) = 2L>t", (0 < a < 2) 



(2) 



then the heat conductivity is given by 



K = — 



-j/VT cx i'' 



(3) 



where the exponent /3 is 



= 2-2/a. 



(4) 



This relation connects heat conduction and diffusion, quantitatively. As 
expected, normal diffusion (a = 1) corresponds to the size-independent 
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{(3 = 0) heat conduction obeying the Fourier law. Moreover, a ballistic motion 
(a = 2) implies that the thermal conductivity is proportional to the system 
size L, a superdiffusive behavior(l < a < 2) implies a divergent thermal con- 
ductivity (/3 > 0), and a subdiffusive behavior (a < 1) implies zero thermal 
conductivity (/3 < 0) in the thermodynamic limit. 

The simple relation (4) is in good agreement with existing data from bil- 
liard gas channels to nonlinear lattices, and even single walled nanotubes[6]. 
This is shown in Fig. 5, where we compare the theoretical prediction (4) with 
existing data in different models. 

We should mention here the one dimensional hard point gas model for 
which anomalous heat conduction has been found by several authors [20, 
21]. However it seems there is no agreement on the numerical value of the 
exponent /3. Indeed in Ref.[20], the value (3 « 0.25 has been found, while the 
value /3 « 0.33 is reported in [21]. 

4 Control of heat flow 

While in the previous sections we have discussed the relation between dynam- 
ical chaos and heat conductivity, in the following we will turn our attention 
to the possibility to control heat flow. Actually a model of thermal rectifier 
has been recently proposed[13] in which the heat can flow preferentially in 
one direction. Although this model is far away from a prototype realization, 
it is based on a mechanism of very general nature and, as such, is suitable of 
improvement and may eventually lead to real applications. This problem is 
discussed in the next section. 

4.1 Thermal diode 

In a recent paper[14], a thermal diode model has been proposed in which, even 
though the underlying physical mechanism is similar to the one in Ref.[13], 
there is a new crucial element which allows to improve the efficiency by more 
than two orders of magnitude. 

The diode model consists of two segments of nonlinear lattices coupled 
together by a harmonic spring with constant strength hint (see Fig. 6). Each 
segment is described by the (dimensionless) Hamiltonian: 



The two ends of the system arc put into contact with thermal baths at 
temperature and Tr for left and right bath, respectively. In fact, Eq. (5) is 
the Hamiltonian of the Prenkel-Kontorova (FK) model which is known to have 
normal heat conduction [23]. For simplicity we set the mass of the particles 
and the lattice constant m = a = 1. Thus the adjustable parameters are 
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Fig. 6. Configuration of tfie tfiermal diode model based on two coupled FK chains. 



(^L, fcini, fcfl, Vl, Vr, Tl, T/f), whcrc the letter L/R indicates the left/right 
segment. In order to reduce the number of adjustable parameters, we set 
Vr = AVl, kR = Afci, Tl = Toil + A), Tr = To(l - A) and, unless otherwise 
stated, we fix = 5, /cl = 1 so that the adjustable parameters are reduced 
to four, (Z\, A, fcini. To)- Notice that when A > the left bath is at higher 
temperature and vice versa when Z\ < 0. 

In Fig. 7 we plot the heat current J versus A for different temperatures 
Tq. It is clearly seen that when A> Q the heat current {,]+) increases with Z\, 
while in the region Z\ < the heat current (J_) is almost zero, i.e. the system 
behaves as a thermal insulator. The results in Fig. 7 show that our model has 
the rectifying effect in a wide range of temperatures. The rectifying efficiency, 
defined as |J+/J_|, could be as high as few hundreds times, depending on 
temperature as well as on other parameters. 




Fig. 7. Heat current J versus the dimensionless temperature difference A for differ- 
ent values of Tq. Here the total number of particles A'' = 100, kint = 0.05, A = 0.2. 
The lines are drawn to guide the eye. 
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Fig. 8. Spectra of the two particles at the interface for different temperatures at 
kint ~ 0. (a) particle at the left side of the interface, (b) particle at the right side 
of the interface. Here A = 0.2, N = 100. 



Rectifying mechanism 

To understand the underlying rectifying mechanism, let's start from the en- 
ergy spectrum of the interface particles. Fig. 8 shows the phonon spectra of 
the left and right interface particles at different temperature when the two 
lattices are decoupled {hnt — 0). 

The match/mismatch of the energy spectra of the two interface particles 
controls the heat current. It is clearly seen from Fig. 8 that, if the left end 
is in contact with the high temperature bath T/,, and the right end with 
the low temperature bath T/j (< 7/,), then the phonon spectra of the two 
particles at interface overlap in a large range of frequencies, thus the heat 
current can easily go through the system from the left end to the right end. 
However, if the left end is at lower temperature Tl and the right end is 
at higher temperature Tr (> Tl), then the phonon spectrum of the right 
interface particle is mainly in the low frequency part, while the left interface 
particle is in the high frequency part. Then there is almost no overlap in 
phonon frequency, and the heat current can hardly go through from right to 
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left, and the system behaves as an insulator. Why the left and right particles 

at the interface have so different phonon spectra? This can be understood 
from the following analysis on temperature dependent phonon spectra, due 
to the nonlinearity. 

(1) Low temperature limit. At low temperature, the particle is confined in 
the valley of the on-site potential. By linearizing the equation of motion one 
can easily obtain the frequency band[14]: 



For the case of Fig.7 with T=0.01 (left) and T=0.002 (right), this corresponds 
to 0.36 < a;/27r < 0.48 for the left particle and to 0.16 < w/27r < 0.21 for the 
right particle. 

As the temperature is increased, the interparticle potential becomes 
more and more important until a critical value Tcr ~ V/{2tt)'^ is reached (we 
take the Boltzman constant equal unity), when the kinetic energy is large 
enough to overcome the on-site potential barrier. At this point low frequency 
appears and this happens at the critical temperatures Tcr = 0.13 for V = 5 
(left), and Tcr = 0.025 for y = 1 (right). This is in quite good agreement 
with the data of Fig. 8. 

(2) High temperature limit. In the high temperature limit the on-site po- 
tential can be neglected, the system is close to two coupled harmonic chains, 
and the phonon band is is given by [24]: 



which gives < u)/2Tr < 0.32 for the left particle and < w/27r < 0.14 for 
the right particle, again in good agreement with Fig. 8. 

In fact, in order to optimize the rectifying effect, one should avoid the 
overlapping of the phonon bands in the low temperature limit (Eq.6) and 
that in the high temperature limit (Eq.7) for each segment of the system. 
According to the above estimates, one should have V > 4fc, which is satisfied 
for the case of Fig. 8. 

We should stress that it is the nonlinearity of the potential that makes 
phonon spectra temperature dependent, and thus the rectifying effect possi- 
ble. Therefore it is reasonable to expect the rectifying effect to be present, in 
different degrees, in typical nonlinear lattices. 

4.2 Negative Differential Thermal Conductance 

Apart from the "one-way heat flow" , the negative differential thermal resis- 
tance phenomenon observed in a certain temperature intervals in the thermal 
diode is of particular interest. As illustrated in Fig.7 for A < —0.2, a smaller 
temperature difference (A), can induce a larger heat current since, due to 
nonlinearity, it can result in a better match in phonon bands. 

The same phenomenon is shown for different inter-face coupling constants 
hint in Fig. 9(a), and different system size in Fig. 9(b). 



< u> < VV + 4k. 



(6) 



0<oj < 2Vk, 



(7) 
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Fig. 9. (a) Heat current versus temperature Tl (at fixed Tr = 0.2) for diflferent 
coupling constants, ki„t, with lattice size A'^ — 50. The system parameters are: 
Vl = 5, Vr — 1, kL = 1, kR = 0.2. (b) Same as (a) but for different system size N. 
kint = 0.05. Notice that when Tl < 0.1 the heat current increases with decreasing 
the external temperature difference. 



4.3 Thermal Transistor 

The phenomenon of negative differential thermal resistance allows us to pro- 
pose a "thermal transistor" [15]. The configuration of the thermal transistor 
is shown in Fig. 10(a). It consists of three segments, D, S and G. The names 
D, S and G follow the ones used in a MOSFET (Metal-Oxide-Semiconductor 
Field-Effect-Transistor) that is the most important device for very large scale 
integrated chips such as microprocessors and semiconductor memories. Seg- 
ment D (from D to O) has a negative differential thermal resistance in a 
certain temperature regime while segment S is a normal heat conductor, i.e., 
heat current inside this segment is positively dependent on temperature dif- 
ference. Segment G is the control segment, which is connected to the junction 
particle between segments S and D. Temperature Tq will be used to control 
temperature To (at the junction O) so as to control the heat current from 
D to S. In analogy to the MOSFET, in which the electronic current in the 
electrode G is very small, we require here that the heat current Jq through 
segment G to be as small as possible, (otherwise it is hard to set Tq to a re- 
quired value in experiment). Moreover the heat resistance of segment G must 
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be small enough in order to well control the temperature To by changing Tq 

so that To-Tg- 

Notice that, in typical situations, the differential heat resistance, Rg = 

I Wr' ) ill segment 5*, and Ru = — ( lip- ) in segment D, 

\oio jrp^ =const \'^^" / Td =const 

are both positive and therefore there exists only one value of Tq for which 
Js = Jd so that Jg = 0, Since Jg = Jd + Jg, the "current amplification 



factor", a 



dJp 

dJc 



< 1, namely in order to make a change AJu, 



Rs+Rd 

the control heat bath has to provide a larger AJq- This means that the 
"transistor" can never work! 

The key point of our transistor model is the "negative differential heat 
resistance" as we observed in the diode model[14]. It provides the possibility 
that when To changes both Js and Jd change simultaneously in the same 
way. Therefore Js = Jd (or Jg w Jd) can be achieved for several different 
values of Tg or even in a wide region of To as shown in Figs. 10 and 11. In this 
situation heat switch and heat modulator/amplifier are possible. In the ideal, 
limiting case of Rs = —Rd which, in principle, can be obtained by adjusting 
parameters, the transistor works perfectly. 



Thermal Switch 



We first demonstrate the "switch" function of our transistor, namely we show 
that the system can act like a good heat conductor or an insulator depending 
on the control temperature. This is illustrated in Fig. 10(b), where we plot 
Jg, Js, and Jd versus Tq- When Tg increases from 0.03 to 0.135, both Jd and 
Js increase. In particular, at three points: Tq ~ 0.04, 0.09 and 0.135, Jd = Js 
thus Jg is exactly zero. These three points correspond to "ofF' , "semi-on" and 
"on" states, at which Jd is 2.4 x lO"'', 1.2 x lO"'^ and 2.3 x lO"'^, respectively. 
The ratio of the heat current at the "on" state and that at the "off' state is 
about 100, hence our model displays one important function - switch - just 
like the function of a MOSFET used in a digital circuit. 



Thermal modulator /amplifier 

As demonstrated above, the heat current from D to S can be switched be- 
tween different values. However, in many cases, like in an analog circuit, we 
need to continuously adjust the current Js and/or Jd in a wide range by 
adjusting the control temperature Tg- In Fig. 11 we demonstrate this "mod- 
ulator/amplifier" function of our transistor. The basic mechanism of such 
"modulator/ amplifier" is the same as that of the "switch" but we consider 
here different parameter values. It is seen that in the temperature interval 
Tg € (0.05,0.135), the heat current through the segment G remains very 
small ((— 10~^ ~ 10~^), within the shadow strip in Fig. 10, while the heat 
currents Jg and Jq continuously increase from 5 x 10~^ to 2 x 10~^. 
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Fig. 10. (a) Configuration of the thermal transistor, (b) Heat current versus the 
control temperature To- Parameters are: Td ~ 0.2, ,Vd ~ 1.0, fco ~ 0.2, ki„t ~ 
0.05; Ts = 0.04, Vs = 5,ks = 0.2, Va = 5,kG = 1, hntc = 1- Notice that both Js 
and Jd inrease when the temperature Tq is increased. 




0.04 0.08 -p 0.12 0.16 



Fig. 11. Heat current versus the control temperature To- Here: Td = 0.2, Vd = 
1.0, ko = 0.2, Ts = 0.04, Vs = 5,ks = 0.2, h„t = 0.05, Va = 5,kG = 1, h„tG = 0.1. 
The shadow region is the range of variation of Jg in the temperature interval 
Tg G (0.05,0.135). 
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5 Conclusions and discussions 

In this paper, we have given a brief summary of our recent work on heat 
conduction in one dimensional systems. We have shown that strong chaos is 
sufficient but not strictly necessary for the validity of the Fourier heat law. 
Indeed linear mixing can be sufficient to induce a diffusive process which 
ensures normal heat conductivity. 

For systems with total momentum conservation one typically finds anoma- 
lous conductivity, namely the thermal conductivity is divergent with the sys- 
tem size. Anomalous conductivity has been connected with anomalous diffu- 
sion via the very simple formula (4). 

Finally we have shown the possibility to build a thermal diode which 
exhibits a very significant rectifying effect in a very wide range of system pa- 
rameters. Moreover, based on the phenomenon of negative differential thermal 
resistance observed in the thermal diode, we have built a theoretical model 
for a thermal transistor. The model displays two basic functions of a transis- 
tor: switch and modulator/amplifier. Although at present it is just a model 
we believe that, sooner or later, it can be realized in a nanoscale system ex- 
periment. After all the Prenkel-Kontorova model used in our simulation is a 
very popular model in condensed matter physics[25]. 



Acknowledgement 

BL is supported in part by Faculty Research Grant of National University of 
Singapore and the Temasek Young Investigator Award of DSTA Singapore 
under Project Agreement POD0410553. GC is partially supported by EU 
Contract No. HPRN-CT-2000-0156 (QTRANS) and by MURST (Prin 2003, 
Ordine e caos nei sistemi estesi non lineari: strutture, stocasticita' debole e 
trasporto anomalo). 

References 

1. R. E. Pcicrls, Theoretical Physics in the Twentieth Century, edited by M. Fiera 
and V. F. Weisskopf (Wiley, New York, 1961). 

2. F. Bonetto et al, in "Mathematical Physics 2000," A. Fokas et aJ. (eds) (Im- 
perial Collogo Press, London, 2000) (pp. 128-150); S. Lepri et al., Phys. Rep. 
377, 1 (2003). 

3. H. Kaburaki and M. Machida, Phys. Lett. A181, 85 (1993);S. Lepri et al. 
Europhys. Lett. 43, 271 (1998);S. Lepri, Phys. Rev. E 58, 7165 (1998); S. 
Lepri et al, Phys. Rev. Lett. 78, 1896 (1997); B. Hu, B. Li, and H. Zhao, Phys. 
Rev. E 61, 3828 (2000); A. Pcrcvcrzcv, Phys. Rev. E, 68, 056124 (2003). 

4. J.-S Wang and B. Li, Phys. Rev. Lett.92, 074302 (2004); Phys. Rev. E, 70, 
021204 (2004). 



Title Suppressed Due to Excessive Length 



15 



5. B. Li and J. Wang, Phys. Rev. Lett. 91, 044301 (2003); 92, 089402 (2004). 

6. B. Li, J. Wang, L Wang, and G Zhang, CHAOS, xxx(2005); cond-mat/0410355. 

7. G. Casati, J. Ford, F Vivaldi, and W. M. Visschcr, Phys. Rev. Lett, 52, 1861 
(1984). 

8. D. Alonso, R. Artuso, G. Casati, and I. Guarneri, Phys. Rev. Lett 82, 
1859(1999). 

9. B. Li, L Wang, and B Hu, Phys. Rev. Lett. 88, 223901 (2002) 

10. D. Alonso et ah, Phys. Rev. E 66, 066131 (2002). 

11. B. Li, G. Casati, and J. Wang, Phys. Rev. E 67, 021204 (2003) 

12. B. Li, G. Casati, J. Wang, and T. Prosen, Phys. Rev. Lett. 92, 254301 (2004) 

13. M. Terraneo, M. Peyrard, and G. Casati, Phys. Rev. Lett 88, 094302 (2002). 

14. B. Li, L. Wang, and G. Casati, Phys. Rev. Lett. 93, 184301 (2004). 

15. B. Li, L Wang, and G. Casati, cond-mat/0410173 

16. B. V. Chirikov, Phys. Rep. 52, 263, 335(1979). 

17. G. Casati and T. Prosen, Phys.Rev. Lett. 85, 4261 (2000). 

18. G. Casati, and T. Prosen, Phys. Rev. Lett. 83, 4729 (1999) 

19. A. Dhar and D. Dhar, Phys. Rev. Lett. 82, 480 (1999). 

20. A. Dhar, Phys. Rev. Lett., 86, 3554 (2001); G. Casati and T. Prosen, Phys. 
Rev. E 67, 015203(R) (2003). 

21. P. Grassberger, W. Nadler, and L. Yang, ibid. 89, 180601 (2002); H. Li and H. 
Zhao, ibid 89, 079401 (2002). 

22. J. Bardeen and W. H. Brattain, Phys. Rev. 74, 230 (1948). 

23. B. Hu, B. Li, and H. Zhao, Phys. Rev. E 57, 2992 (1998). 

24. C. Kittel, Introduction to Solid State Physics, 7'th edition, John Wiley & Sons, 
Inc. (New York), 1996. 

25. O. M. Braun and Y. S. Kivshar, Phys. Rep. 306, 1 (1998). 



